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Generalizations of the classic Bell inequality to higher dimensional quantum systems known as 
qudits are reputed to exhibit a higher degree of robustness to noise, but such claims are based on 
one particular noise model. We analyze the violation of the Collins-Gisin-Linden-Massar-Popescu 
inequality subject to more realistic noise sources and their scaling with dimension. This analysis 
is inspired by potential Bell inequality experiments with superconducting resonator-based qudits. 
We find that the robustness of the inequality to noise generally decreases with increasing qudit 
dimension. 
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I. INTRODUCTION 

A Bell inequality [T] experiment consists of n parties 
who share an entangled state of n (or more) particles. 
Each party chooses to perform one of m measurements, 
with each measurement producing one of d outcomes. 
By repeating this round multiple times, the probabilities 
for the various joint outcomes can be estimated. A Bell 
inequality for this scenario is a relationship that these 
probabilities must satisfy, if they arise from a local real¬ 
istic model. The general structure of these inequalities 
has been studied intensively since Bell’s original argu¬ 
ment [2]. 

The most famous form of the Bell inequality is the 
Clauser-Horne-Shimony-Holt (CHSH) inequality [3] for 
n = 2 parties (e.g. Alice and Bob), each party perform¬ 
ing one of m = 2 measurement choices, with each mea¬ 
surement registering one of d = 2 outcomes. In terms of 
the joint probabilities, the CHSH inequality reads 

p{Ai = Bi)—p{Ai = H2)+p(A2 = Hi)+p(A2 = B 2 ) < 2, 

( 1 ) 

where the measurement settings are labeled by 1 and 2 
for Alice’s (or Bob’s) choice of measurement, with out¬ 
comes Ai and A 2 for Alice (and Bi and B 2 for Bob), 
and, in a slight abuse of notation, the joint probabil¬ 
ity p{Aq = Bq) indicates the probability that Alice and 
Bob’s measurement outcomes are identical. For certain 
entangled states and measurement choices, quantum me¬ 
chanics predicts, and experiments confirm, a violation of 
this inequality [T]. 

A common experimental procedure, illustrated in Fig. 

replaces the alternative measurement settings by a uni¬ 
tary transformation chosen by the parties and performed 
just before a fixed measurement. This unitary must be 
chosen and performed by one party sufficiently quickly 
to ensure that no information can propagate to the other 
party. If this cannot be done, the experiment is subject to 
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FIG. 1: A general framework for a Bell inequality experiment 
with n = 2 parties. A source of entanglement is shared to 
Alice and Bob, who choose to rotate their part of the joint 
quantum state by unitaries Ha, a and respectively. The 

measurement outcomes Ma and Mb are then recorded, com¬ 
piled, and compared against the inequality. 


the so-called locality loophole, for which the derivation of 
Eq. 0 fails. Another requirement is that the measure¬ 
ments be sufficiently accurate so that the probabilities 
entering the inequality can be reliably estimated. If this 
is not the case, the experiment is subject to the so-called 
detection loophole. Modern experiments with photons 
typically close the locality loophole HE], but are subject 
to the detection loophole (in the guise of the fair-sampling 
assumption, but see recent progress 013). Experiments 
with matter qubits (atoms [8], ions 0 EO], solid-state 
spin qubits [in ED and superconducting qubits [13 E]) 
typically close the detection loophole, but are subject 
to the locality loophole. Achieving a loophole-free Bell 
inequality experiment would be a landmark test of quan¬ 
tum mechanics. 

Towards this goal, and to better understand entangle¬ 
ment and nonlocality in general, there have been many 
studies of generalized Bell inequalites [T]. These include 
the Mermin inequality [15] for multiple qubits (n > 2), 
the Collins-Gisin inequality [T6| for multiple measure¬ 
ments (m > 2), and the Collins-Gisin-Linden-Massar- 
Popescu (CGLMP) inequality [17] for higher-dimensional 
systems known as qudits (d > 2). A recent approach [18] 
to the Collins-Gisin inequality for entangled qudits (with 
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n = 2 and m = d > 2) exhibits the potential to reduce 
the requirements to close the detection loophole. 

In this paper we focus on understanding the CGLMP 
inequality (for n = m = 2 and d > 2), which takes the 
form 

E [l-j^A[V{k)-V{-k-l)]<2, (2) 

/c=0 ^ ^ 

where 

V{k) = P{Ai =Bi^k)p P{Bi =A2Pkpl) 

PP{A2 = B2Pk)p P{B2 = + k). (3) 

Here the joint probabilities are defined for outcomes Aa = 
0,l,...,d — 1, and the addition is performed modulo d. 
These can expressed as 

d-l 

P{Aa = Bb + k) = '^P{Aa = j,Bb = j + k mod d). (4) 

j=0 

We have also studied a closely related inequality proposed 
by Zohren and Gill [19] 

P{A2 < B 2 ) + P{B2 < Ai) + P(Hi < Pi) 

+P(Pi < A 2 ) ^ 1. (5) 

These two inequalities have the remarkable property that 
the violation increases with increasing d. Alternatively, 
if the initial state is subject to depolarizing noise, the 
amount of noise that removes the entanglement, and 
hence the potential to violate the inequality, increases 
with the dimension [20]. In this sense, these higher¬ 
dimensional Bell inequalities exhibit a surprising robust¬ 
ness to noise, and may be useful when exploring advanced 
tests of quantum systems. 

A recent experiment tested the CGMLP inequality us¬ 
ing orbital angular momentum (0AM) states of light [21], 
but did not find the enhanced violation with dimension. 
In fact, they found that the violation ceased for dimen¬ 
sions higher than d = 12. This fact strongly suggests that 
the robustness of the inequality can be compromised by 
the actual noise subject to the system. Previous work 
on this topic has focused only on a simple form of de¬ 
polarizing noise [171120], and its impact on the detector 
efficiency needed to observe a violation of the inequality 
[ 22 ]. To account for the behavior seen in the recent ex¬ 
periment, and to predict the possible violations in other 
experiments, more realistic noise models are necessary. 

In general, noise can affect all of the stages of the 
experiment: entangled-state preparation, state rotation, 
and state measurement. The noise from each stage must 
be analyzed to determine the robustness of the inequality. 
In this paper we examine a general framework for testing 
the GGLMP inequality with qudits, analyze the poten¬ 
tial complexity of each stage of the experiment, and eval¬ 
uate the impact of different types of noise (depolarizing, 
dephasing, and amplitude damping) on the inequality. 



FIG. 2: A quantum circuit for a Bell inequality experiment, 
consisting of entangled-state preparation, single-qudit rota¬ 
tions, and qudit-state measurement. In the first stage, two 
qudits, initially in the |00) state, are acted on by generalized 
Hadamard gates H, a controlled-phase gate (with 0 = 271 jd) 
to prepared a maximally entangled two-qudit state. The 
single-qudit rotations Ua and Ub determine the two measure¬ 
ment bases (depending on the choices of Alice and Bob), while 
the actual measurement is in a fixed qudit basis. 

Our work is inspired by theoretical proposals to use the 
quantum states of superconducting resonators as qudits 
[ 23 II 25 ] ; such systems have the potential for long-distance 
entanglement through microwave photons How¬ 

ever, our analysis is intended to be general enough to 
guide experiments with other potential matter qudit sys¬ 
tems, and our results may have implications for photonic 
qudits as well. Our analysis indicates that the GGLMP 
inequality does not generally exhibit the special robust¬ 
ness to noise claimed in previous work. 

This paper is organized as follows. Section H describes 
a general Bell inequality experiment with ideal qudits, 
analyzing how state preparation, rotation, and measure¬ 
ment can be implemented for a general qudit system. 
Section HI analyzes the potential inequality violation 
subject to depolarizing, amplitude-damping, and dephas¬ 
ing noise on the qudits during the full experimental se¬ 
quence. In Section IV we conclude and discuss outstand¬ 
ing questions. 


II. BELL INEQUALITY EXPERIMENT WITH 
IDEAL QUDITS 

In this section we will describe an analysis of a bipartite 
{n = 2) Bell inequality experiment, focusing on the gen¬ 
eral structure of entangled-state preparation, single-qudit 
rotations, and qudit-state measurement. A quantum cir¬ 
cuit for this process is shown in Fig. [^ We will analyze 
each component of this circuit, drawing inspiration from 
recent theoretical work on the control of superconduct¬ 
ing resonators. However, we expect our results can be 
applied to alternative qudit implementations. 


A. Entangled-State Preparation 

The first important stage of a Bell inequality experi¬ 
ment is the production of entangled states of two systems. 
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As shown in Fig. this can be accomplished by using 
a generalized Hadamard or DFT (for Discrete-Fourier- 
Transform) operation on each qudit and a controlled- 
phase gate between the qudits. 

The generalized Hadamard gate is the discrete Fourier 
transform (DFT) matrix with matrix elements Hj^k = 
where uu = For d = 4, this has the 

matrix form 


H{d = 4) = DFT^ 


1 

71 


/I 

1 

1 

Vi 




(6) 


In general this matrix has d independent elements (and 
total elements), so an analysis of control on qudits 
suggests that implementing a general qudit Hadamard 
will require a number of elementary operations that is a 
polynomial in d [28] . 

The controlled-phase gate is a two-qudit gate that im¬ 
plements the transformation 

cn{e)\j,k) = e-^^^’‘^^\j,k). (7) 


This is a natural generalization of the two-qubit con¬ 
trolled phase gate, and can also be implemented in a 
time polynomial in d [28] . 

The combination of these operations in the preparation 
stage of Fig. [^performs the transformation 


1^) 


{H O I)Cn{2TT/d){H ® H)\0, 0 ) 

. d-1 

il 0 H)Cn{ 2 n/d)- \j^k) 

j,k=0 

. d-1 

j,k=0 


J 


(8) 


This maximally entangled two-qudit state will be used in 
the Bell inequality measurements for Alice and Bob. 

We note that there are alternative approaches to gen¬ 
erating the initially entangled state for superconducting 
resonators, which typically scale linearly in the qudit di¬ 
mension [29] . Qudit operations can also be optimized to 
scale linearly with qudit dimension m- We optimisti¬ 
cally conclude that the state-preparation stage can be 
performed in a timescale that is linear in d. 


B. Single-Qudit Rotations 


rotating their half of the joint-qudit state by one of the 
unitary operators UA,a oi* These have the matrix 

elements 


[UA,a]3,k = ( 9 ) 

and 

[UB,b]3,k = ( 10 ) 

where ai = 0, 0^2 = 1/2, j3i = 1/4, and P 2 = —1/4. 
These operations involve the discrete Fourier transform 
(or its inverse), along with diagonal phases. These again 
can be performed in a time polynomial in d, which we 
again optimistically take as linear in the qudit dimension 

m- 


C. state Measurement 

In order to close both the locality and detection loop¬ 
holes, one wants a fast and efficient measurement of the 
qudit states. For some qudit implementations, such as 
hyperfine states of atoms or multilevel superconducting 
devices, one can implement a direct d outcome measure¬ 
ment. For qubit-resonator systems in trapped ion, cavity, 
or circuit-QED systems, one often implements an indi¬ 
rect measurement coupling the resonator to one or more 
qubits. We consider a state measurement approach that 
maps the resonator state of d = 2^^ dimensions onto n 
qubits, inspired by the discussion in Chapter 6 of m- 
After this mapping, one can measure the state bit by bit. 
This allows a single-shot measurement of the qudit state. 

The measurement scheme can be represented by a 
quantum circuit, illustrated in Fig. composed of single¬ 
qubit Hadamard gates, two-qubit controlled-phase gates, 
and a special qubit-resonator controlled-phase gate. This 
last gate acts on a qubit-resonator state as follows: 

|^)qubit C) |^)res ^ ^ ^|^)qubit C) |^)res- 

This gate can be implemented using the natural evolu¬ 
tion of a qubit-resonator system in the dispersive regime 
[24] . resonant qubit-resonator logic operations [32], or 
other methods [25l [33]. In brief, the measurement 
circuit performs the mapping by sequentially correlat¬ 
ing the resonator state \y) with a register of qubits 
({xi,X 2 ,... ,Xn}). Using a bit string {yi ,... ,^^1 to la¬ 
bel the resonator state \y), with the (big-endian) binary 
representation 

d = ^ 7 -y2‘^^ ^ + • • • + dn-l2 + dn? (12) 

the circuit performs the mapping 

|00 • • • 0)qubits C) |d)res ^ |dld2 ' ' ’ dn)qubits ^ |d)res- (13) 


After Alice and Bob have chosen their measurement 
basis (a or 6), they will adjust their measurement by 


Then, measurement of the qubits will provide a bit-by-bit 
measurement of the resonator state \y). 
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FIG. 3: A measurement circuit to map the quantum state of a resonator (of dimension d — 2'^) onto the state of n qubits. This 
circuit is composed of 2n single-qubit Hadamard operations, n qubit-resonator controlled-phase gates (with phases indicated in 
the circular gates on the resonator), and n(n+l)/2 qubit-qubit controlled-phase gates (with phases indicated in the rectangular 
gates on the qubits). 


We now analyze the measurement circuit of Fig. [^in 
more detail. The qubits are initially prepared in an equal 
superposition state (by Hadamard gates): 

^ H Cy\xx---Xn)®\y) (14) 

{xx,...,Xn},y 

The first qubit-resonator controlled-phase gate produces 
the state 

^ H Cy{-lf^y\xi---Xn)®\y). (15) 

{xi,...,Xn},y 

Using the binary representation, we have 


The next qubit-resonator phase gate produces the con¬ 
trolled phase 



(- 1 ) 


Xn-iyn-l 


(i) 


Xn— 1 yn 


( 20 ) 


where we have used an argument similar to Eq. (16). 
The final phase can be eliminated by a qubit-qubit phase 
gate, while a subsequent Hadamard on qubit n — 1 yields 


^ Cy{-ir-^y-^+--^^\x^---Xr,_^zyr,)®\y). 

{xi,...,Xn-i,z},y 

( 21 ) 

Performing the sums over Xn-i and z, using the corre¬ 
sponding form of Eq. ( p!8| ) , the second stage yields 


(_2)31nl/ 


^_ '^^Xjiy\‘2 ^_ 'Y^^'ny2‘^ . . . ^_ 'Y^Xjiyji 

[-ly^yry^ (16) 


A subsequent Hadamard on qubit n yields 


1 ^ 2 ) 2(«-2)/2 E 

{xi,...,Xn-2},y 

After k stages we will find 


Xn-2yn-iyn) ® \y)■ 

( 22 ) 


^ Cy{-ir-y-+--^\x,---Xn-iz)®\y). 

{xi,...,Xn,z},y 

(17) 

Performing the sums over Xn and 2 ), using the fact that 


^ ^ (18) 

Xri =0 

we find 

1^1) = 2(n-l)/2 E Cy\xx---Xn-iyn)®\y)- 

{xx,...,Xn-i},y 


2 {n-k)l 2 ^ 

^ ^ Cyl^i • • • ^72—/c^n—/c-|-l ’ ' ’ Vn) G) 1^) 

{xx,...,Xn-k},y 

(23) 

The next qubit-resonator gate produces the controlled 
phase 

^gi 7 r/ 2 ^^ ^ _ -^\^Xn-kyn-k ^ 

/ . ,^\ Xn-kyn-k + 1 / . ,^k\^'ri-kyn 

(e^^/2\ .. /gt 7 r /2 \ 


( 19 ) 


( 24 ) 
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Eliminating all of the trailing controlled phases by qubit- preserving map 
qubit phase gates, performing a Hadamard gate on qubit 

n — k, and evaluating the summation over x^-k yields Pn+i = E EmPnEly^^ (27) 

m 

l^/c+i) = 2 (n-/c-i )/2 ^ where the error operators Em depend on the type of 

noise, to be defined below, and are parametrized by a 
/ ^ Cy\xi • • • Xn-k-iVn-k ''' Vn) ^ Iv) • single number p. This map is iterated N times, where 
{xi,...,xn-k-i},y N is proportional to the number of fundamental steps 

(25) during the preparation and rotation stages of the Bell 
circuit of Fig. and p is proportional to the fidelity of 


Repeating until /c = n, we finally obtain 

l^'n) = ■■■2/") ® l2/)- (26) 

y 

A measurement of the n qubits will produce the outcome 
(^ 1 , ^ 2 , • • •, ^n) with probability \cy\^. This circuit uses 
n qubit-resonator gates, 2n Hadamard gates, and n{n — 
l)/2 qubit-qubit controlled-phase gates. Thus, we see 
that this circuit requires a timescale of order (log 2 d)^. 


III. BELL INEQUALITY EXPERIMENT WITH 
NOISY QUDITS 

The preceding analysis shows that, in general, one can 
expect that the state preparation, rotation, and measure¬ 
ment stages of a Bell inequality experiment will each scale 
with the qudit dimension. If these stages are composed of 
fundamental elements, such as quantum gates, that are 
each subject to noise, then the resulting experiment will 
be subject to noise that scales with the qudit dimension. 
How that scaling affects the inequality is the subject of 
this section. 

While we will analyze idealized models of noisy qudits, 
our results can be understood in physical terms. Each 
fundamental quantum gate will take some time, during 
which the quantum state can be subject to fluctuating 
fields, lose energy to the environment, or lose quantum 
information in some other fashion. Each of these noise 
processes will affect the quantum state of the system, 
in terms of its density matrix. If the total time for the 
experiment scales with the qudit dimension, the result¬ 
ing density matrix will be subject to a correspondingly 
increased amount of noise. 

A detailed analysis of this process would require model¬ 
ing the qudit s’ Hamiltonian and its coupling to external 
fields and the environment. To understand the scaling 
with dimension, however, we can simplify our analysis 
to consider a single parameter figure-of-merit for each 
step of the experiment, such as the gate fidelity, and look 
at three types of noise: depolarizing noise, amplitude 
damping, and dephasing noise [34]. Using these models, 
we allow the strength of the noise to vary with the qudit 
dimension in the following way. 

We take the initial density matrix po = |^)(^|5 with 
|T) given by Eq. (§, and iterate the appropriate trace- 


each step. The final density matrix, after N applications 
of the noise map and subject to the measurement choices 
a and 6, is given by 

Pa,6 = {UA,a ^ UB,b)pN{UA,a ^ - (28) 

The probabilities of the measurement outcomes are then 
P{Aa =j,Bb = k) = {j,k\pa,b\j^k). (29) 

As discussed above, the number of steps required to 
prepare the initial state and rotate the measurement is 
taken to be linear in d, while the single-shot measurement 
circuit scales with (logd)^. For the dimensions d < 16 
considered below, there is little difference between linear 
and logarithmic scaling. Thus, to understand how scaling 
affects the Bell parameter we set N = d. For com¬ 
parison, we will also consider A = 1, similar to previous 
work. 


A. Depolarizing Noise 

Depolarizing noise simulates the interaction of a sys¬ 
tem with a high temperature environment. With each 
iteration, the system has probability 1 — p of becoming 
depolarized, i.e. replaced by the completely mixed state, 
so that 

Pn+l=PPn + {i-p)-^, (30) 

while the system has probability p of being unaffected. 
Note that for this map, smaller values p indicate stronger 
noise (this convention is chosen to match [II]). As de¬ 
scribed above, we allow this map to be repeated for a 
variable number of times before calculating the Bell pa¬ 
rameter. 

The value of the generalized Bell parameter Id{p) is 
shown in Fig. as a function of d and for several values 
of the probability p. From these results, we can make 
several observations. First, for p = 1, we see that the Bell 
parameter increases with dimension, as found previously. 
In fact, we find 

= I) » 2.97(1-5^). ( 31 ) 

Second, for p < 1 and A = 1, the Bell parameter again 
increases with dimension, as seen in Fig. |^a). However, 
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Depolarizing Noise: No Scaling Depolarizing Noise: Linear Scaling 




FIG. 4: Generalized Bell parameter Id as a function of qudit dimension d for two types of depolarizing noise and various noise 
strengths, (a) Single application of the depolarizing operator, with probabilities of 1, 0.99, 0.95, and 0.9, from top to bottom, 
(b) Repeated application of the depolarizing operator with the number of applications is linear with the qudit dimension, with 
probabilities of 1, 0.99, 0.95, and 0.9, from top to bottom. 


the overall violation does decrease with decreasing p. In 
fact, it can be shown that, for this type of depolarizing 
noise, Id{p) = pld{^) [II]- This behavior underlies the 
reputed robustness to noise. 

By constrast, if we apply the map N = d times before 
calculating Id [Fig. [^b)], we find that the violation de¬ 
cays with dimension. This can be understood by noting 
that repeated iteration of Eq. (30) has the solution 


The specific model of amplitude damping is given by 
the map 


1 

Pn+l = E {Ee(E>Em)pn{Ee(S)Emy, (33) 

i,m=0 

where the single-qudit amplitude-damping operators are 
given by 


Pd=p‘^Po + (1 


(32) 


SO that, in the presence of depolarizing noise scaling with 
d, the Bell parameter behaves as p^Id{p = !)• Thus, the 
small increase in the Bell parameter (for p = 1) of Eq. 
(31) is quickly reduced by the exponentially decreasing 
factor p^. 

This may have some bearing on the experimental re¬ 
sults discussed above [21] . If one of the stages of the ex¬ 
periment has a fidelity that decays with the dimension of 
the entangled state, as evidenced in previous experiments 
[35] . one can easily reproduce the observed decay of the 
inequality violation. Eor example, setting p = 0.998 and 
N = d produces results in rough agreement with Eig. 3 
of El]. 


B. Amplitude Damping Noise 

Amplitude damping describes the effects of energy dis¬ 
sipation on the quantum system. We use a simplified am¬ 
plitude damping model, chosen to represent the damping 
of a quantum resonator over a time At with decay time 
T. In this model, the singly excited state |1) decays with 
probability p = while state \j) decays with prob¬ 

ability jp [^ . 


d-l 

= ( 34 ) 

i=o 

and 

d-l 

Ei = ^yr^ij-i)(ii. (35) 

i=i 

We note that while this model intrinsically scales with 
dimension, we continue to allow the number of iterations 
of the map to scale with dimension as well. It is also 
relevant to observe that this map is an approximation 
to the actual decay process, in that single application of 
this map with p = 0 corresponds not to complete loss of 
energy, but to the subtraction of a single photon. 

Using this model of amplitude damping, the resulting 
Bell parameter Id as a function of dimension d is shown 
in Eig. for several values of the probability p. Eor 
A = 1, the violation decays slowly with d, but remains a 
violation for up to d = 16 for p = 0.9. Already, however, 
we see a significant difference from the behavior of Id 
under depolarizing noise. 

Eor N = d iterations, the violation decays much faster 
with d. However, the decay is surprisingly less than the 
depolarizing map for the same value of p. We will soon 
return to this issue, but for certain values of d and choices 
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of the probability, the two are similar. For example, to 
model the 0AM experiment [21] , one could use amplitude 
damping with p = 0.992 and N = d. 

The results above have used a discrete model of energy 
loss, in which a single quantum is removed from the sys¬ 
tem with some probability. Of course, the actual physics 
involves a continuous loss of energy. We have compared 
our discrete amplitude damping model to a more contin¬ 
uous model, by breaking up the decay into a number of 
intervals with the same overall decay factor p. As one 
might expect, the difference between the two models is 
negligible for p near 1, in which the probability of losing 
two quanta can be neglected. 


C. Dephasing Noise 

Dephasing describes loss of quantum information with¬ 
out loss of energy: rather than changing the amplitudes 
of the states as a function of time, the energy eigenstates 
of a system accrue random phases with some probability. 
This is modeled by 


d-l 

Pn+l = PPn + il-p) y] {{j,k\Pn\j,k))\j,k){j,k\. (36) 
j,k=0 

After each iteration, the off-diagonal elements of pn are 
reduced by a factor p. 

We note that the off-diagonal elements of the final den¬ 
sity matrix in this case are identical to those found for 
the depolarizing map of Eq. (30). Remarkably, when 
measuring the density matrix using the DFT operators, 
the joint probabilities of Eq. (29) do not depend on the 
diagonal elements. We thus arrive at the interesting re¬ 
sult the effect of dephasing noise is identical to that of 
the depolarizing noise considered previously. 

To verify that the probabilities do not depend on the 
diagonal elements, we consider an arbitrary “diagonal” 
density matrix of the form 


Pdiag — ^ ^ b '5 (i: ^1: 

j,k=l 


(37) 


where 


j,k=l 


(38) 


For this density matrix we use Eq. (29) to calculate the 
joint probablitiies 

77T/) = (8) Ub, bPdiag (UA^a ^ 

= E c^AmAAj)?\{rn\UBAk)? 


j,k=l 

d 


— J_ — J_ 


j,k=l 


(39) 


Since these probabilities are independent of the diag¬ 
onal elements of any density matrix will contribute the 
constant value of l/d^ to the probability P{i^m). Of 
course, the off-diagonal elements will also contribute to 
the probability. However, any two density matrices with 
the same off-diagonal elements will lead to the same prob¬ 
abilities, confirming the result claimed above. 

We can also use this observation regarding the proba¬ 
bilities to help understand the difference between ampli¬ 
tude damping and depolarizing / dephasing noise. For 
a single qubit, amplitude damping reduces the diagonal 
elements by p, but the off-diagonal elements by y/p > p. 
Thus, the effect of amplitude damping on the probabil¬ 
ities will be less than that for depolarizing / dephasing 
noise, for the same value of p. That this holds true for 
dimensions d > 2 is still somewhat surprising. One pos¬ 
sible explanation is that coherence between neighboring 
qudit states is largely maintained, much like the decay of 
a coherent oscillator state m, but we have not explored 
this conjecture further. 


D. Thresholds for Inequality Violation 

A convenient way to summarize the results obtained 
above is to consider the minimum probability Pmin for 
which the CGLMP inequality can be violated. This prob¬ 
ability satisfies 


^d{Pmin) — 2 , ( 40 ) 

and depends on the different models of noise. Recall 
that a smaller probability p indicates a greater amount of 
noise. For p < Pmin, the inequality is not violated, so we 
call Pmin the threshold probability. Previous work found 
that Prnin dccreascs with dimensions, with the conclusion 
that higher dimensional Bell inequalities can be violated 
even if the system is subject to greater amounts of noise. 

We have calculated the minimum probabilities for the 
various noise models described above and their depen¬ 
dence on dimension, as shown in Fig. We find that 
only for depolarizing / dephasing noise with N = 1 does 
the threshold increase with dimension. For all other 
cases—depolarizing / dephasing noise with N = d, and 
amplitude damping with N = 1 and N = d —the min¬ 
imum probability increases with dimension. For these 
noise models, a violation of the CGLMP inequality with 
higher-dimensional systems requires higher-fidelity oper¬ 
ations. 


E. Alternative Procedures 

We have studied a number of alternative procedures 
of the Bell inequality test. First, we have looked at the 
inequality proposed by Zohren and Gill [T9| described in 
the introduction, and find that the threshold probabilities 
are identical to those shown in Fig. [^ Second, we have 





Amplitude Damping: No Scaling 



(b) 

Amplitude Damping: Linear Scaling 



FIG. 5: Generalized Bell parameter Id as a function of qudit dimension d for two types of amplitude damping and various 
noise strengths, (a) Single application of the amplitude damping operators, with probabilities of 1, 0.99, 0.95, and 0.9, from 
top to bottom, (b) Repeated application of the amplitude damping operators, with the number of applications is linear with 
the qudit dimension, with probabilities of 1, 0.99, 0.95, and 0.9, from top to bottom. 


CGLMP Violation Thresholds 



FIG. 6: The minimum probability Pmin foi' which the GGLMP 
inequality is violated, as a function of qudit dimension d and 
for the noise models (from top to bottom) depolarizing / de¬ 
phasing noise with N — d^ amplitude damping with N — d^ 
amplitude damping with V = 1, and depolarizing / dephas¬ 
ing noise with iV = 1. Systems with noise p < Pmin will not 
violate the inequality. 


looked at the alternative entangled states 


where A/* is a normalization factor; these states can 
achieve a higher degree of violation than the maximally 
entangled states [37], and approximate those states with 
maximal violation. While the thresholds are slightly dif¬ 
ferent for these states, they exhibit the same general be¬ 
havior seen in Fig. [^ 

Finally, we have considered another variation of the 


entangled state 


^ d—l 

= ( 42 ) 

^ ^ j=0 

in which Bob’s qudit states have been reversed (Bob’s 
unitaries must also be “reversed”). These states 

are easier to produce in superconducting circuits [29|, 
and one might think they would be less sensitive to am¬ 
plitude damping. However, the thresholds are again only 
slightly different for this alternative procedure. In short, 
the behavior seen in Fig. appears to be generic for 
these inequalities and noise models. 


IV. CONCLUSION 

We have studied the effects of noise on Bell inequal¬ 
ity experiments using the higher-dimensional inequalites 
proposed by Collins, Gisin, Linden, Massar, and Popescu 
m and Zohren and Gill [19] and variations thereof. By 
modeling the required operations needed in the prepa¬ 
ration, rotation, and single-shot measurement stages 
of these experiments, we have analyzed how the num¬ 
ber of operations scale with the qudit dimension. For 
most types of noisy operations, namely amplitude damp¬ 
ing, depolarizing, and dephasing noise, we find that 
the higher-dimensional inequalities require increasingly 
higher-fidelity operations. This conclusion runs counter 
to previous work, which considered depolarizing noise 
with a fixed noise probability, independent of qudit di¬ 
mension. However, this conclusion is in agreement with 
a very recent analysis of higher-dimensional Bell inequal¬ 
ities with random coherent errors [38] . 

We conclude with a few observations. Our results show 
that these higher-dimensional inequalities are not more 
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robust against noise when taking into account the full ex¬ 
perimental procedure. Thus, these inequalities do not ap¬ 
pear to offer a quick route towards closing the detection 
loophole with photons. However, there are other multi¬ 
setting {m > 2) inequalities with qudits [18] that allow 
for reduced detection efficiencies for both atom-photon 
and photon-photon entanglement, and in the presence of 
depolarizing noise. Our results motivate continued anal¬ 
ysis of these and other inequalities (including those with 
n > 2) in the presence of realistic noise appropriate for 
matter and photonic qudits. The continued acquisition 


of evidence for nonlocality from a wide variety of physi¬ 
cal systems remains an intriguing and important goal of 
modern quantum physics. 
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